Mechanical properties of carbynes investigated by ab initio total-energy calculations 
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As sp carbon chains (carbynes) are relatively rigid molecular objects, can we exploit them as 
construction elements in nanomechanics? To answer this question, we investigate their remarkable 
mechanical properties by ab-initio total-energy simulations. In particular, we evaluate their linear 
response to small longitudinal and bending deformations and their failure limits for longitudinal 
compression and elongation. 

PACS numbers: 62.25.-g, 81.07.Gf, 62.20.de, 46.70.Dc 
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I. INTRODUCTION 

The rich chemistry of carbon is due to the capa- 
bility of its electronic configuration to adjust to dif- 
ferent bonding situation. Carbon atoms realize three 
main hybridization schemes of the valence orbitals: sp 3 , 
sp 2 , and sp hybridization, sp 3 bonding tends to form 
three-dimensional networks, as in diamond and the re- 
lated amorphous structures. sp 2 -hybridized orbitals lead 
to two-dimensional networks, at the basis of graphene, 
graphite, fullerenes, nanotubes, ribbon structures, and 
other forms of wide current interest for their unique me- 
chanical and electronic properties, sp-bybridized car- 
bon forms linear structures (carbynes), which, compared 
to their chemically more stable sp 3 (alkanes) and sp 2 
(alkencs) hydrogenated counterparts, tend to be rigid, 
thus potentially appealing as backbones for molecular 
nanotechnology. 

Carbynes of increasing length and varied termina- 
tions are being synthesized^— and characterized £ At 
the same time, the fabrication of carbon atomic chains 
from stretched nanotubes or graphene was achieved by 
controlled electron irradiation in transmission electron 
microscopes^— The remarkable robustness of these car- 
bon chains under irradiation combined with the ease of 
electron-beam fabrication at the nm scale can provide 
a route to the synthesis of actual nano-devices based 
on carbynes. Recent investigations have elucidated the 
prominent role of molecular oxygen as a primary source of 
chemical degradation of sp carbon^S. In a pure-carbon 
environment carbynes are far more stable^— . 

Effects of longitudinal strainingi£r— and lateral 
bending^ were investigated very recently by ab-initio 
methods. The outcome of those studies is that carbyne 
chains are at one time extremely stiff against longitudinal 
straining and very soft against a bending deformation, to 
the extent that even extreme bending affects only mod- 
erately their bonding properties^ In the present paper 
we focus on the mechanical properties of carbynes in the 
perspective of exploiting them as construction materials 
for nano- engineering. We use ab-initio simulations^ 3 -^ 




FIG. 1: (Color online) A beam model used to describe the 
mechanical properties of a sp-carbon chain subject to strain- 
ing. The scheme represents the relaxed configuration of the 
beam subjected to four lateral constraints, which induce two 
bending moments on the beam central section, which, as a 
result, bends as the arc of a circle. We identify the end-to- 
end distance l s of the deformed region; the curvilinear beam 
length /; the maximum deflection d from straight geometry; 
the radius of curvature r; and the deflection angle cj>. 



to predict the stiffness of the carbynes for both longitu- 
dinal strain and for bending, and their ultimate tensile 
strength. The main results of these simulations is that a 
sp carbon chain (and it matters little whether it is of cu- 
mulenic - all double bonds, or polyynic - with alternating 
single-triple bonds) exhibits not only an exceptional me- 
chanical stiffness against longitudinal deformations, but 
also a small but nonzero rigidity against bending. Ac- 
cordingly, a C n chain resembles a thin beam character- 
ized by a finite rigidity against buckling, rather than a 
string. We find that, e.g., a single-carbyne Cs "pillar" can 
withstand a modest but non-negligible axial compression 
force ~ 0.7 nN before buckling. 



II. THE MODEL 

In simulations, one can tune the carbyne length by fix- 
ing the position of the end carbon atoms of a selected 
chain section, and leave all other atoms free to relax to 
their equilibrium position, compatible with the imposed 
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strain. As well know to civil engineers, as soon as the 
compressive stress exceeds a critical value, the straight 
configuration of a beam subject to a purely axial strain 
e < turns unstable and "buckles" spontaneously to a 
curved arc. In the simulations of highly strained con- 
figurations, it is straightforward to address either the 
metastable straight configuration, or the stable buckled 
configuration investigated in Ref. l20l 

To characterize the mechanical properties of a carbyne 
segment, and in particular to evaluate the buckling tran- 
sition, one needs to identify the correct linear-response 
parameteres in a simple model where the effects of differ- 
ent kinds of deformations can be decoupled and studied 
separately. The simplest continuum model describes a 
carbyne as a strained thin beam of a homogeneous mate- 
rial, thereby ignoring its discrete atomic structure. Three 
basic deformations can be applied to the beam: longitudi- 
nal traction or compression; lateral bending, as sketched 
in Fig. [TJ torsion around its axis. Correspondingly, the 
elastic deformation energy can be decomposed as: 
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where E tcns is the tensile energy increase due to longi- 
tudinal elongation/shortening, -Ebend is the energy due 
to the pure flexion of the beam, and -Btors is the energy 
increase due to the twisting of the beam. Remarkably, 
even this third term is important for the essentially ID 
cumuleneSiH In the present work we will however mostly 
ignore torsion, and concentrate on the first two terms. 

In the linear-response (small-strain) regime, we write 
the tensile energy in terms of the strained length I relative 
to the fully relaxed length Iq as 



-Etcns - "ft' - ^o) 2 



(2) 



where x/^o is the elastic stiffness of the beam, and \ is the 
length-independent beam elastic tension: in the contin- 
uum mechanics of 3D solids x = £ • A, the product of the 
material's elasticity modulus times the beam transverse 
area. 

Figure Q] indicates the relevant geometrical quantities 
characterizing the pure-bending deformations of a car- 
byne of fully relaxed length Iq. In a circular-arc geometry 
ensuing from the application of a small bending moment 
M at the beam ends, the bending energy can be derived 
by the linear relation of the deflection angle 4> with the 
applied bending moment: 
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Here g is the length-independent bending stiffness of the 
beam: in 3D continuum mechanichs g = E ■ I, the prod- 
uct of the elasticity modulus times the second moment 
of area of the beam cross section. (The dimension of g 
is force x area.) For a circular deformation, the bending 
energy -Ebend is therefore 



^bend = \;M<t> = 9 
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(4) 



It is convenient to express all quantities in terms of 
clearly defined geometrical quantities, e.g. the distance l s 
between the end points of the chain and the the maximum 
deflection d from straight geometry. We have: 
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With the implicit substitution of the expressions ([5]) and 
(O, the elastic energy of the beam stretched and/or bent 
in a circular arc is therefore 
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III. CALCULATIONS 

We determine the linear-response parameters \ an< i 9 
by means of simulations based on the density functional 
theory (DFT) in the local density approximation (LDA). 
The time-honored LDA is one in many functionals be- 
ing used for current DFT studies of molecular systems: 
other functionals often improve one or another of the sys- 
tematic defects of LDA (underestimation of the energy 
gap, small overbinding and overestimation of the vibra- 
tional frequencies), but to date no functional is univer- 
sally accepted to provide systematically better accuracy 
than LDA for all properties of arbitrary systems. For a 
covalent system of s and p electrons as the one studied 
here, LDA is appropriate, and we expect our results to 
change by a few percent at most if the calculations were 
repeated using some other popular functional^ - — 

We compute the total adiabatic energies by means of 
the plane-waves DFT code Quantum Espresso . 25 39 We 
consider the two different limiting structures of carbynes: 
polyynes (or a— carbynes) with alternating single-triple 
bonds (. . . — C=C— C=C— . . . ), contrasted to cumulcncs 
or (3— carbynes, with nearly-equal-length (double) bonds 
(. . . =C=C=. . . ). A typical polyynic chain is obtained 
when each end carbon atom forms a single bond with a 
ligand (e.g. a hydrogen atom, as in diacetylene) , while 
a typical cumulene is obtained when each end carbon 
forms a double bond, for example to a CH 2 group, as in 
ethylene. 

We focus on a Cs chain segment, and compare its me- 
chanical response in its HC2-C8-C2H polyynic and H 2 C- 
C8-CH2 cumulenic realizations. For the polyyne, we se- 
lect hydrogenacetylide terminations rather than a single 
hydrogen because the HC2- group allows us to impose 
a bending moment at the ends of the Cs chain as in 
Fig. [TJ by displacing a C atom along the chain continua- 
tion rather, than a chemically different H atom, thereby 
probing the intrinsic properties of carbyne, rather than 
those of a specifically-terminated compound. We also 
consider C n polyyne of different lengths, to check for 
size effects, and a 90° end-twisted Cs cumulene, which 
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FIG. 2: (Color online) The DFT-LDA tensile excess energy 
computed for the strained straight HC2-C8-C2H polyyne (bul- 
lets), as a function of e — (l — lo)/la- Solid line: a quadratic fit 
of Eq. ([2}, restricted to the 7 data points within ±0.8% strain, 
where deviations from the linear-response regime are negligi- 
ble. Crosses: the total excess energy without the constraint 
that the chain should remain straight: buckling starts to build 
up for strain exceeding e cr it- The very similar data and fit for 
the cumulene H2C-C8-CH2 is not reported for clarity's sake. 



we must treat within the local spin-density approxima- 
tion (LSDA) since twisting induces a peculiar total spin- 
1 electronic state, with a twofold-degenerate level occu- 
pied by two parallel-spin electrons^ We do not attempt 
any serious size scaling, for two main reasons: (i) very 
long carbynes are mechanically unstable (thus unsuitable 
for nanomechanical applications) anyway, and (ii) the 
long-range interaction effects demonstrated recently^ 
would make a proper ab-initio size-scaling determination 
of the mechanical (in particular bending) properties of 
long carbynes prohibitively expensive. In the present 
work we do not deal with the exotic odd-n C„ chains, 
which are chemically less stable and more difficult to 
synthesize i 26 i 27 

To evaluate Xi starting from the fully relaxed carbyne 
we stretch or compress a straight Cs chain section by fix- 
ing the positions of the end carbon atoms, and letting all 
other atoms free to relax to their resulting equilibrium 
positions under stress. We extract an estimate of x by 
fitting Eq. @ to the resulting values of the excess en- 
ergy as a function of the chain elongation, as reported in 
Fig.H 

To evaluate the bending stiffness g, we force the chain 
to bend by imposing a small lateral displacement of 
the two atoms adjacent to the Cs section of the chain, 
while keeping the first and last atom in the chain sec- 
tion bounded to remain laterally undisplaced, following 
the scheme of Figs. [1] and [3] (inset). We then let all other 
degrees of freedom (including the axial position of the lat- 
erally bounded atoms) of the chain relax. The Cs section 
relaxes to form a circular arc, which we fits** to extract 
the geometrical quantities indicated in Fig. [T] Note that 
this procedure differs substantially from the fixed-end ap- 
proach of Ref. |2(X where longitudinal and bending strains 
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FIG. 3: (Color online) The bending energy Eh end — E mGc h 
Etens, see Eq. (0, as a function of the squared bending angle 
of a -Cs- chains. Here -Emcch is the excess energy in the bent 
configuration induced by the 4-point fixed lateral displace- 
ment (at most 30 pm) sketched in Fig. [1] iJtcns is the value 
obtained by Eq. @, using \ as obtained in the straight-chain 
calculation of Fig. [2] Symbols are the DFT values, lines are 
linear fits, whose slopes, according to Eq. Q, represent the 
values of g. Polyynic chain: circles and dashed line; in-plane 
bent cumulenic chain: squares and dotted line; the out-of- 
plane bent cumulenic chain results in very similar energies 
and is not shown for clarity. Inset: an example of fully relaxed 
in-plane bent configuration of the cumulene H2C-C8-CH2. 
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FIG. 4: (Color online) The chain-length dependency of the 
computed bending stiffness g, as reported in Table [I] 



are applied at the same time, and deviations extend well 
outside the linear-response region. In the language of 
that work, the maximum arc-chord ratio in our calcu- 
lations is l/l s = 1.0064. We repeat this procedure for 
several values of the fixed lateral displacements. With 
the computed DFT energies and structures, we extract 
the values of g by fitting the angular dependency of -Eb en d 
in Eq. (gj), see Fig.O 
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(cumulene) 


(polyyne) 
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(90°-twisted) 






la (C„) [pm] 


377 


886 






889 


1395 


1903 


£ult 


16.2% 


16.3% 


17.9% 


16.8% 


16.4% 


16.8% 


F u it [nN] 


12.7 


12.3 


14.7 


12.7 


12.1 


12.1 


X [nN] 


150 


156 


151 


148 


149 


148 








out-of-plane 


in-plane 








g [nN(nm) 2 ] 


7.91 1(T 2 


5.69 1(T 2 


5.68 10" 2 


5.95 10~ 2 


5.62 10~ 2 


5.11 10 -2 


5.03 10" 2 


Ecrit 


-3.7% 


-0.458% 


-0.470% 


-0.493% 


-0.473% 


-0.174% 


-0.093% 


Fcrit [nN] 


5.49 


0.715 


0.710 


0.744 


0.702 


0.259 


0.137 



TABLE I: The computed mechanical properties of selected carbynes. The elastic tension \ and bending stiffness g define the 
elastic response, Eq. ([7]). The ultimate strain e u i t and tension F u i t define the maximum stretching that the beam can stand 
before it breaks. For a Cs chain section with both ends pinned (hinged, free to rotate), F cr it measures the minimum compressive 
force (with a corresponding strain E cr it = (?Jj rlt — Iq)/Io relative to the reported equilibrium length lo) that must be applied to 
induce buckling. 



The resulting elasticity parameters x an( i 9 are c °l~ 
lected in Table Q] Despite the different bonding configu- 
ration of polyyne and cumulene, the overall longitudinal 
elasticity x values are remarkably similar. The value of 
X is compatible with a sound velocity v s = (x/a*) 1 ^ 2 — 
31.5 km/s for long-wavelength longitudinal vibrations (p, 
is the linear mass density). In view of the long-range 
interactions demonstrated recently^, the weak depen- 
dence of x on t ne chain length is remarkable. The bend- 
ing stiffness g is comparatively more sensitive to the 
chain length, see also Fig. 2J If this velocity is used to 
evaluate the highest-frequency Brillouin-zone-boundary 
phonon in a simplest mass-spring chain, one obtains a 
frequency h*BZb — v s /(7rdcc) — 79 THz, i.e. a wave num- 
ber ~ 2600 cm -1 , expectedly slightly above the region of 
the observed high-frequency Raman and infrared modes 
emerging as the characteristic spectroscopic signatures of 
the carbynes ^MtM 

The cumulene is anisotropic: one can make it bend ei- 
ther within the molecular plane, as sketched in the inset 
of Fig. [31 or perpendicularly to it, resulting in differently 
oriented double bonds^S being perturbed. We find that 
bending is significantly softer for the out-of-plane distor- 
tion than for the in-plane one, with a ~ 5% anisotropy. 

The proposed cosine expression for the strain energy 
given in Eq. (9) of Ref. [20 is compatible with our Eq. (gj) 
for small curvature (although a small amount of lon- 
gitudinal strain energy is neglected there). The rela- 
tion g ~ dccFd obtained based on the average bond 
length dec = 129 pm and orientation-dependent bond- 
ing strength Fd — 2.54 eV as obtained in Ref. yields 
g = 5.26 10 -2 nN/(nm) 2 , in good agreement with our 
result, especially with the C12 chain value. This agree- 
ment is even more surprising in view of the extreme de- 
formations introduced in Ref. [20I and indicates that the 
linear-response coefficient of Eq. Q extends well outside 
its validity region, provided that the cosine espression 



Eq. (9) of that work is adopted. 

According to elementary continuum mechanics^ the 
critical axial load that a beam can sustain in its straight 
configuration before the onset of the buckling instability 
is given by Euler's formula 

^crit = I , (8) 
corresponding to a critical strain 

^ = ^.-^-.3 * (9) 
'0 X X'o 

for a critical sinusoidal lateral deformation profile. By 
substituting the computed values for x an d <?, we find 
that the Cg carbyne sustains a critical strain of nearly 
—0.5%, under the action of a critical force of approxi- 
mately 0.7 nN, before buckling, as reported in Table |U 
The corresponding critical strain and force for cumu- 
lene (in the softer out-of-plane direction) are essentially 
equal. The computed values are relevant when both ends 
of Cs are pinned (hinged, free to rotate), while if both 
ends were frozen, the critical force and strain would be 
four times larger. Likewise, following Eqs. © and ©, 
-F cr it decays rapidly with size. For example, for Ci6, 
whose length is approximately 2.15 times that of Cs, 
-Fcrit = 0.14 nN, i.e. approximately 20% of the Cs value 
only. Direct calculation done for the C12 polyyne in HC2- 
C12-C2H agrees with this scaling, with a small deviation 
due mostly to a smaller value of g, relative to the Cs 
value (while x is practically coinciding, see Table [I]). 

IV. DISCUSSION AND CONCLUSION 

We find it rather surprising that, despite the presence 
of in principle softer single bonds, the polyynic chain is 
essentially as hard to compress as double-bonds-based 
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cumulene. The same essentially equivalent stiffness of 
polyync and cumulene is found against bending defor- 
mations, but here the result is less unexpected. Even 
a maximally twisted cumulene in a high-spin electronic 
state exhibits very similar mechanical properties, see Ta- 
ble m 

We check the obtained buckling critical point for possi- 
ble effects of chain discreteness or anharmonicity, by run- 
ning actual simulations. We compare pairs of DFT simu- 
lations at fixed strain: one bound to the linear configura- 
tion as in Fig. [5J and one starting off in a slightly curved 
geometry. In the curved simulations, below the buck- 
ling instability, the relaxing chain goes back to straight, 
and recovers the same energy and outward force that 
the chain produces on the pinned end atoms as in the 
straight calculation. In contrast, above the buckling in- 
stability, a curved shape stabilizes, with a net decrease 
in total energy and in the force acting on the pinned 
end atoms. Crosses in Fig. [5] report the fixed-end ex- 
cess energy of the polyynic chain when allowed to relax 
in a buckled geometry. We find that the actual buck- 
ling instability occurs very close to the linear-response 
continuum-model value e cr ;t ofEq. Relaxations near 
instability are rather delicate, since equilibrium is almost 
indifferent, providing a manifold of almost-equivalent ge- 
ometrical configurations which give a hard time to the 
optimization algorithm. 

The present work demonstrates that it is possible to 
go a long way in describing an immaterialy thin ob- 
ject such as a monoatomic carbon chain with the me- 
chanics of bulk construction elements, not unlike it was 
done for 2D graphene in Refs. [32ll33l Eventually, car- 
byne chains turn out rigid enough that the 886 pm-long 
free-standing Cs polyyne can sustain a compressive strain 
of nearly 0.5% before buckling. On the tensile side, we 
extend our calculations well outside the linear-response 
region, to estimate the ultimate tensile strength of car- 



bynes. The ultimate tension F u it is computed as the max- 
imum force that the carbyne section produces in sustain- 
ing an externally imposed longitudinal strain e u u > 0. 
The result in the 10 nN region agrees with previous 
determination s 17 ' 34 (but disagrees significantly with the 
Tersoff-Brenner-model determination of Ref . 1351) . and in- 
dicates that a single molecular chain could be used as a 
rope to lift a mass as heavy as one microgram without 
breaking! If one could pack carbynes with a lateral den- 
sity of one per A = 0.2 x 0.2 nm 2 cross-section, one would 
obtain a material characterized by a remarkable ultimate 
tensile strength of the order 300 GPa, comparable to that 
of carbon nanotubes^ 

Despite the remarkable mechanical properties of car- 
bynes, their chemical reactivity makes them unsuitable 
to many applications in real world mechanical situations, 
where carbon fiber or carbon nanotubes provide superior 
stability with comparable mechanical properties. How- 
ever in clean well-isolated nano-engineered devices one 
could envisage that the usage of C„ chains as structural 
or elastic elements may lead to consistent advantages over 
traditional solutions. For example, the soft bending de- 
gree of freedom could be exploited for the construction of 
sensitive accelerometers. Indeed, according to our evalu- 
ation of g, a 0.1 /xg mass hanging at the end of a 37 nm- 
long C30 chain (with the other end bonded to a fixed 
substrate) would deflect laterally by as much as 1 nm 
under an acceleration of 0.01 ms -2 , i.e. one thousandth 
of the Earth gravitational field. 
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